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Chapter  VII 

and  Spherical  Premixed  Planes 


1.  Cylindrical  Flames 

Although  it  has  not  been  studied  to  the  same  extent  as  the  plane  premixed 
flame,  the  cylindrical  flame  is  in  principle  almost  as  easy  to  produce.  The 

reacting  mixture  is  supplied  through  the  surface  of  a  circular  cylinder  and  is 

* 

induced  to  flow  radially  by  means  of  sufficiently  close  end  plates.  The  flame 
then  forms  a  coaxial  cylinder  and  can  be  observed  through  the  end  plates,  which 
should  be  transparent  and  good  thermal  insulators. 

Analytically  the  cylindrical  flame  stands  between  the  plane  and  spherical 
flames.  The  structure  of  its  reaction  zone  is  the  same  as  that  of  the  plane 
flame,  with  temperature  constant  beyond;  so  that  there  is  no  curvature  effect 
as  for  the  spherical  flame.  On  the  other  hand,  like  the  spherical  flame  it 
does  not  exhibit  the  cold-boundary  difficulty:  the  mixture  must  be  introduced 
at  a  finite  radius,  which  can  however  be  so  small  that  a  line  source  is 
effectively  formed.  Ironically  enough,  in  their  attempt  to  treat  curved 
flames  Spalding  &  Jain  (1959)  use  plane-flame  analysis  on  the  spherical  flame, 
where  it  is  never  valid,  and  neglect  the  cylindrical  flame,  where  it  is  always 
valid.  fhe  results  are  qualitatively  correct,  however 

(Ludford,  1976). 

The  object  of  Secs.  2  &  3  is  to  show  how  cylindrical  geometry  modifies 
a  premixed  flame  (Ludford  1980).  For  simplicity  we  shall  consider  a  single 
reactant  (monopropellant)  undergoing  a  first-order  reaction  and  take  lC  =  1. 


Spherical  premixed  flames,  which  have  many  similarities  to  spherical  diffusion 
flames,  will  be  treated  in  Secs.  U— 6 . 
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2 ,  Planar  Character 

We  shall  only  consider  the  basic  features  of  cylindrical  flames;  the 
discussion  can  easily  be  extended  to  other  questions  treated  for  plane  flames 
in  Chs.  II  &  III.  When  making  comparisons  with  the  latter  bear  in  mind  that  the 
non-diraensionalization  is  here  based  on  the  radius,  a,  of  the  supply  cylinder. 

If  v  is  the  radial  velocity,  the  equation  of  continuity  admits  rpv  being 
constant.  When  M  is  the  mass  flux  at  the  supply  cylinder  we  may  therefore 
write 


(1) 


rpv  =  M  or  pv  =  M/r, 


the  latter  giving  the  mass  flux  at  every  other  radius  r.  Once  the  temperature 
T  has  been  determined,  the  density  p  =  1/T  and  v  =  KT/r  follow  immediately, 
while  the  small  variations  in  pressure  about  its  constant  level  can  be  calculated 
from  the  momentum  equation  of  the  mixture.  There  remain  then  the  energy  and 
reactant  species  equations 


(2) 


jC  (Y)  =  -£( T)  =  dy  exp(-e/T), 


where 


(3) 


/  =Il(r  L.)  JL 

'*’•  r  drv  dr  r  dr 


in  the  cylindrical  geometry. 

First  note  that  the  Shvab-Zeldovich  variable  Y  =  Y  +  T  satisfies 


00 


( Y)  =  0  for  1  <  r  <  <». 


Since  the  only  solutions  which  remain  bounded  at  infinity  are  constants,  we 
have 
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(5)  Y  +  T  =  T  =  Y  ‘  +  T 

'  00  s  s 

for  a  reaction  that  goes  to  completion,  where  s  denotes  conditions  at  the 
supply  r  =  1.  The  fact  that  the  Y,  T-relation  is  identical  that  for 

plane  flames  (Sec.  II. 2)  is  responsible  for  the  structures  of  the  reaction  zones 
being  the  same. 

The  asymptotic  analysis  of  the  equations  (2)  proceeds  as  for  a  plane  flame. 
The* temperature  beyond  the  flame  sheet  must  be  constant,  i.e. 

(6)  Y  =  0,  T  =  T^  for  r  >  r#, 

while  up  to  the  flame  sheet  the  reaction  is  frozen,  i.e.  /(Y)  =  /(T)  =  0, 
so  that 

(7)  Y  =  Y  +  L(1  -  rM),  T  =  T  +  L(rM  -  1)  for  1  <  r  <  r#. 

s  s 

» 

Here  L  =  M-1T' ,  with  T'  the  temperature  gradient  at  the  supply,  is  the  heat 
s  s 

conducted  back  into  the  supply  per  unit  mass  of  mixture.  These  two  pairs  of 
formulas  (6),  (7)  give  the  same  values  at 

(8)  r*  =  [1  +  (Tm  -  Ts)/L]1/M; 

as  expected,  the  stand-off  distance  for  a  plane  f] |$ne  is  recovered  as  the  radius 
a  of  the  supply  cylinder  tends  to  infinity  when  due  attention  is  paid  to  a 
mass-flux  unit  proportional  to  l/a.  Consistency  requires  r#  >  1,  i.e. 

(9)  I.  >  0, 

which  means  the  supply  must  be  a  conductive  heat  sink. 

As  before  the  interior  of  the  flame  sheet  is  investigated  with  the  expansions 


(10) 


Y  =  6y(f,)  +  0(5),  T  =  +  6t(,£)  +  o(6) 


vhere 


(11)  C  =  (r  -  rJ/5  and  6  =  T2/0. 

"  09 

The  structure  is  thereby  found  to  satisfy  the  equation 

(12)  d2t/d£2  =  -D  y  e\  vhere  y  +  t  =  0, 
and  the  boundary  conditions 

(13)  t  =  KJgC/r#  +  o(l)  as  ?  -v  -to  ,  t  =  o(l)  as  5  4  +», 

which  come  from  matching  with  the  expressions  (6)  and  (7)  outside.  Here 
?  -0/T 

(lit)  D  =  6  e  "D  while  J  =  Y  -  Y’/M 

s  s  s 

is  the  reactant  flux  fraction  Y  -  rY'/M  at  the  supply  (usually  l).  Exactly 

the  same  problem  is  obtained  for  a  plane  flame  except  that  the  coefficient 

MJ  /r#  in  the  condition  (13a)  is  replaced  by  (which  equals  MJ  there), 

s  s 

Noting  that  here  J#  =  J  /r#  (because  the  total  reactant  flux  2-nrVJ  is 

s 

conserved  for  frozen  chemistry)  shows  that  the  cylindrical  reaction  zone  is 
locally  plane.  For  spherical  flames  the  y,t-relation  (12b),  which  derives 
from  the  Y,T-relation  (5)  is  changed;  so  that  the  reaction  zone  is  quite 
different  from  its  plane  counterpart. 

From  the  solution  (ll.2l)  of  the  corresponding  plane-flame  problem  we  deduce 

(15)  D  =  (Js/r#)2!.^/2 
so  that 

p  p  6/T  h  9  9 

(16)  D  =  ( J  0  e  “/2T44)ir, 

which  is  the  required  M,D-relation.  (The  presence  of  r#  in  this  formula  means 
that  it  depends  on  all  supply  parameters,  unlike  the  plane  flame.)  For  fixed 


ML,  i.e.  heat  conducted  back  to  the  supply,  it  has  the  general  shape  of  the 
parabola  (II.22)  obtained  for  the  plane  flame,  because  the  factor 

(17)  1/r*  =  [1  +  (Tm  -  T  )/L]"2^M 

so  5 

varies  only  between  exp[-2(T  -  T  )/ML]  and  1  as  M  increases  from  0  to  ». 

co  5 

The  parabola  is  useful  for  determining  the  speed  (i.e.  M  )  with  which  a  plane  flame 

propagates  into  fresh  mixture  at  given  pressui’e  (i.e.  D) ,  but  there  is  no 

equivalent  use  here.  On  the  other  hand,  for  the  set-up  envisaged  in  the  opening 

paragraph  both  M  and  D  are  prescribed  (along  with  T  and  J  ),  and  the 

s  s 

formula  determines  the  final  temperature  T  (note  L=J+T-T). 

»  s  s 

3«  Near-Surface  and  Surface  Flames.  Remote  Flames. 

For  the  above  solution  to  be  valid  the  parameter  values  must  be  such  that 
r#,  as  given  by  the  formula  (8),  lies  between  1  and  ®  .  As  for  the  plane  flame 
(Sec.  II. 5)  three  limiting  cases  arise,  two  of  whiqh  are  essentially  the  same 
as  there.  The  third  leads  to  an  interesting  new  phenomenon. 

When  M  and  D  become  large,  with  all  other  parameters  held  fixed,  the 
flame  approaches  the  surface.  An  intense  convective-diffusive  zone  of  thickness 
0(M  )  forms  near  the  surface,  bounded  by  a  reaction  zone  whose  thickness  is 

ote^M”1). 

By  contrast  a  true  surface  flame  can  be  produced  for  any  M  by  adjusting 
the  pressure  so  as  to  make  T^  -*■  Tg.  Similarly  remote  flames  can  be  produced 
by  making  T^  ■*  Jg  +  Tg  (i.e.  L  -»■  0).  At  both  extremes  the  preceding  analysis 
becomes  invalid:  either  the  boundary  intrudes  into  the  reaction  zone  and  there 
is  no  frozen  region  between  them  or  the  isothermal  limit  of  (7)  is  not  uniformly 
valid  in  the  unbounded  frozen  region.  In  either  case  the  asymptotics  must  be 


reworked . 


The  analysis  of  the  surface  flame  is  identical  that  that  for  the  plane  case 
provided  x  is  changed  to  r-1.  We  conclude  that  D  will 
change  from  J2^/2 ,  the  value  (15)  when  r#  =  1,  to  *  as  the  temperature 

S  4 

difference  T  -T  ,  measured  on  the  6-scale,  decreases  from  »  to  0, 

»  s 

By  contrast,  the  remote  flame  cannot  be  treated  as  in  the  plane  case  since 

the  asymptotic  analysis  breaks  down  earlier,  in  fact  as  soon  as  L  becomes  0(d). 

The  difference  lies  in  the  reactant  flux  M  J  /r„  at  the  flame,  which  now  tends 

s  * 

to  zero  like  6^^  as  r*  ->  »;  the  condition  (13a)  loses  its  effectiveness 
unless  a  different  scale  is  used  for  the  structure.  Setting 

(18)  r  =  r*  +  d1"1^ 


gives  the  new  condition 


(19)  t  =  MJ  (£/J  r'  V  o(l)  as  g  -*•  -co. 

s  s 


where  L  =  6£;  the  corr-~  ending  change 
(20)  6  =  i2(1-1/H>  e’8/T”D 


must  also  be  made  to  keep  the  structure  equation  balanced. 

We  therefore  end  with  the  same  problem,  except  that  r*  is  replaced  by 

(j  /£)^^j  so  that  D  is  given  by  the  formula  ( 15 )  with  the  same  replacement  and 
s 

(21)  J2(l-1/M>  e2(l-1/M)  ■ 

S  oo 

Thus  D  varies  from  0  to  »  on  the  scale  e 'v  ”  e  as  j,  increases 


from  0  to  «  ,  the  upper  end  of  the  range  corresponding  to  0  on  the  orevious 

2  e/T« 
scale  o  e 


The  most  interesting  feature  is  the  spreading  of  the  zone  in  which  there  is 
chemical  activity,  as  M  decreases.  The  transformation  (l8)  implies  that  its 
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thickness  is  o(9^k*^) ,  so  that  for  M  <_  1  it  is  no  longer  a  sheet;  indeed 
for  M  <  1  it  has  infinite  extent  [remember  r#  =  0(0“^)  is  larger  still]. 

Such  a  phenomenon  should  be  easily  observable. 

We  have  ssen  that,  except  when  remote,  the  cylindrical  flame  is  locally 
plane,  unlike  the  spherical  flame.  These  results  stem  from  the  diffusion- 
convection  operator 

* 

(22)  /  =  .  M  f-] 

**  a  L3r  3r  3rJ 
r 

governing  the  reactionless  field  behind  the  flame.  Here  a  =  0  (plane,  vhen 
r  =  x),  1  (cylindrical)  or  2  (spherical)  and 

(23)  £  (T)  =  0 
has  the  general  solution 

(24)  T  =  A  +  B  e1*  (a  =  0),  A  +  B  rM  (a  =  l),  A  +  B  e"M^r  (a  =  2) 

in  the  three  cases.  Boundedness  of  T  makes  B  =  0  in  the  first  two  cases 

but  not  in  the  third,  where  T  becomes  an  assignable  parameter  in  addition 

00 

to  any  others.  It  is  this  difference  between  the  convection-diffusion,  process 
in  plane  and  cylindrical  geometries  on  the  one  hand  and  spherical  geometry  on  the 
other  which  accounts  for  the  similarities  and  differences  of  the  corresponding 
flames. 

4.  Spherical  Flames.  Dar.kohler-Nur.ber  Asymptotics. 

If  general  values  of  and  had  been  retained  in  Ch.  VI,  most  of 
the  remaining  results  in  the  present  chapter  could  now  be  oltc.ir.ed  by  setting 
V1  =  0*  v2  =  1  al  =  °»  a2  =  Kornsandia  &  Luaford  19&0).  The 

exceptions  concern  strong  burning  with  M  -»■  »,  However,  the  analysis  need  only  be 


sketched,  omitted  details  being  similar  to  those  in  Ch.  VI. 


We  first  turn  to  the  limits  D  0,  ",  which  had  already  been  considered 
"by  Fendell  (1969).  Frozen  combustion  is  described  by  the  same  formulas  as  for 
the  diffusion  flame,  provided  the  new  definition 

-6/T 

(31)  e  =  De 

is  made;  in  particular  the  burning  rate  (VI.ll)  holds.  Equilibrium  limits 
(p  *  »)  require 

(32)  Y  -  0 

outside  the  reaction  zone,  a  state  that  cannot  hold  next  to  the  source: 
the  zone  must  therefore  lie  on  the  surface  of  the  source.  There  are  two 
possibilities,  according  as  M  is  finite  or  not;  the  former  is  the  counterpart 
of  the  surface  diffusion  flame  and  the  latter  corresponds  to  the  Burke-Schumann 
flame',  as  we  shall  now  see. 

If  general  stoichiometric  coefficients  had  been^ retained,  the  formulas 
(VI.  1^,  15)  would  have  had  -Z  /2a„  in  place  of  Z  .  Then  M  ■»  »  and 
r#  1  as  oip  ->•  0:  the  Burke-Schumann  solution  tends  to  a  surface  flame  with 
indefinitely  large  burning  rate  as  the  oxidant  becomes  inactive.  The  limit 
behavior  cannot  be  deduced  from  generalizations  of  the  expansions  (VI.  17,18,19) 
but  must  be  determined  ab  initio. 

For  that  purpose,  note  that  the  Shvab-Zeldovich  relation  (28)  gives 

(33)  T  =  T  for  r  >  1 

a 

with  an  adjustment  to  T^  at  distances  0(m).  The  temperature  must  therefore 
vary  to  leading  order  in  the  reaction  zone;  as  also  follows  from  the  boundary 
condition  M  ^T*  =  L,  which  in  addition  points  to 
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as  the  appropriate  coordinate  there.  We  conclude  that  the  reaction  zone  ui 
governed  by 

(35)  (d2/dx2  -  d/dx)T  =  DM-2Yexp'(-0/T-) ,  Y  +  T  =  T&, 

(36)  T  -  T  ,  dT/dx  =  L  for  x  =  0  and  T  T  as  x  -*•  », 

S  8. 

i.e.  the  plane  problem  for  a  vaporizing  liquid  considered  in  Sec.  11,3.  It 
follows  that 

(37)  DM"2  =  a(t  ,L,e) , 

s 

where  A  is  the  eigenvalue  introduced  there,  determines  how  M  tends  to 
infinity  with  D. 

The  asymptotic  form  (II.22)  of  A  as  6  ■>  «  was  the  goal  of  Sec.  II. k,  its 
determination  for  finite  0  being  a  numerical  question  on  which  a  great  deal 
of  effort  has  been  spent  in  the  past.  For  the  diffusion  flame  the  Burke- 
Schumann  solution  can  be  written  in  simple  analytical  terms  for  any  0  ,  but  here 
we  end  with  the  corresponding  plane  problem.  Even  the  asymptotic  behavior  of 
M,  the  counterpart  of  the  formula  (VI.l4),  is  a  numerical  question  .  One 
simple  feature  does  remain,  however:  the  plane  solution  only  exists  for 

(38)  L  <  1 

[a  limiting  form  of  the  condition  (VI.2o)],  as  was  shown  in  Sec.  II. b. 

We  turn  now  to  the  equilibrium  limit  which  is  the  counterpart  of  the  surface 
diffusion  flame;  for  convenience  we  shall  refer  to  it  as  the  Buckmaster  limit 
also.  M  remains  finite,  the  temperature  is  constant  to  leading  order  in 
the  reaction  zone  and  the  Shvab-Zeldovich  relation  (28)  gives 

M  =  M  =  ?,n[l  +  (T  -  T  )/(L-l)  ] 

S  00  S 


(39) 
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in  place  of  the  result  (VI.21).  The  expansions  (VI. 22 >23)  are  replaced  by 

(UO)  M  =  M  [1  +  D*1//2/{l~L)  +  oCD^1^2)] 

s 

and 

(41)  T  =  T.  +  d;1/2Ms[1  -  e"5  +  (L-lhl  +  o(d;1/2) 
if  the  new  definition 

-0/T 

(42)  D#  =  De 

is  made.  The  changes  amount  to  2  being  replaced  by  1  to  account  for  the 

different  (non-dimensional)  heat  release  at  the  flame.  Restrictions  on  the 

parameters  come  from  M  ,  which  must  be  positive  if  the  (small)  mass  fraction 

s 

of  reactant  is  to  be  positive: 

(43)  T  >  T  according  as  L  ^  1. 

The  conditions  for  the  existence  of  the  three  limits  determine  three  regions 

in  the  parameter  plane  (Fig.  l)  in  each  of  which  two  of  the  limits  hold.  In  the 

fourth  region  there  is  no  solution  for  any  values  of  D  and  0:  the  Shvab-Zeldovich 

relation  (28)  would  require  Y  to  be  negative  there. 

s 

5.  Ignition  and  Extinction. 

We  turn  now  to  activation-energy  asymptotics,  starting  with  the  nearly 
adiabatic  flame  for  £.r  1,  i.e.  the  analog  of  Sec.  VI. 3.  Again,  different 
values  of  k  in  the  defining  condition  (VlJ$  will  be  generated  by  varying 
T  keeping  T  and  L  fixed  with 

S  co  , 

(44)  L  <  l. 


»■--* ....  j.  .... 


so  that  elimination  of  C  yields  the  basic  equation 


bMdM/dx  =  k  - 


-M/r*  .  0  0  «  0/T  i, 

5  r1  T2  An(M2©2  e  ?2D£r#  T4) 


where  r#(M)  is  the  function  (48).  The  steady-state  response  for  k  =  0 
has  a  simple  algebraic  form  in  which  D  is  proportional  to  (M  +  £n  L)4/M^, 
a  result  essentially  found  by  Fendell  (1972).  Otherwise  the  above  analysis 
has  not  been  published  before. 

Steady-state  responses  are  shown  in  Fig.  2.  They  differ  from  those  for  the 
diffusion  flame  (Fig.  VI. 4)  because  of  the  factor  Z *  there,  which  vanishes 
at  the  Burke-Schumann  value  (VI .6k)-,  here  M  tends  to  infinity  (like  D^2) 
on  all  curves.  For  large  values  of  k  the  turning  points  are 

(52)  K  -  £n(l/L)  +  ljT2/k,  D  -  U/e^lVe^  72(£n  L)2^k4 


(53)  M  ~  fin(k/LT2)  +  An  £n(k/LT2)  -  An  2, 

CO  CO 

0 J  ip 

(54)  D  -  (e2  e  (M2  +  4  An  L  M  +  6  £n2L). 


gf*  *  } 
m,  1 


Conditions  for  instability,  predicted  by  the  full  equation  (5l)  are  the  same 
as  for  the  fuel  drop:  for  <£  <  1  the  middle  branch  of  an  S-response  is 
unstable;  for,,^  >  1  the  other  two  branches,  as  well  as  the  whole  of  a  monotonic 
curve,  are  unstable. 

These  results  and  those  in  the  previous  section  suggest  that  an  S-response 
will  be  obtained  in  the  triangular  region 


T  <  T  <  T 
s  eo  a 


sgy.  i 

ft.  | 


1 

W  \ 


of  the  parameter  plane  (Fig.  l).  The  first  inequality  ensures  that  a  weak- 
burning  limit  exists,  while  the  second  comes  from  the  existence  of  such 
near-adiabatic  responses  for  k  sufficiently  large.  As  for  the  diffusion 
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flame,  the  lower  branch  (including  the  ignition  point)  will  be  obtained  for 
values  M  within  0(9  of  its  frozen  value.  A  similar  analysis  cannot 
be  made  of  the  upper  branch,  however,  because  the  Eurke- Schumann  limit  is 
now  infinite.  The  treatment  required  has  been  given  by  Kapila,  Ludford  & 
Buckmaster  ( 197  5 ) - 

Take  X-  1  again.  At  any  point  on  the  lower  branch  the  temperature 
is  still  given  by  the  formulas  ( VI. 66,67).  Moreover,  at  infinity  the  equation 
(Vl.TO)  and  boundary  conditions  (VI.7l)  still  hold,  provided  only  that 

2  “0/T» 

(56)  D  -  D6~2  e 

V 

now.  The  lower  branch  is  threfore  given  by  Fig.  VI. 6  again,  which  may  be  used 
to  read  off  the  ignition  point. 

Although  M  goes  to  infinity  in  the  equilibrium  limit  the  combustion  field 
for  strong  burning,  i.e.  M  large  but  finite,  has  the  same  structure  (35,36). 
In  the  limit  6  ->  «  it  shows  a  flame  sheet  located  within  o(pf1)  of  the 
surface,  with  equlibrium  beyond.  (Since  the  decomposition  zone  has  thickness 
0(0  1)  on  the  x-scale,  it  stands  clear  of  the  surface  however  large  M  i~.) 

We  therefore  have  the  plane  problem  of  Sec.  II. I*  except  that  the  temperature 
is  not  constant  behind  the  flame  but  varies  according  to  equation  (28)  with 
Y  =  0.  As  a  consequence  the  flame  temperature  is  perturbed  away  from  T0 : 

(57)  T*  »  T  +  (T  -  T  )e“M/L. 

*  a  co  a 

The  result  (II.22),  with  J  =  1  and  T^  replaced  by  this  T#,  gives  the 
response  curve 

(58)  D  =  M2  exp(9e  ^) 


to  leading  order,  where 
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(59)  D  =  [2Tl'e“2exp(-0/T  )]D  and  0  =  [(T  -T  )/LT2]e. 

O.  Q.  a  oo  a 


The  response  has  a  turning  point  at  K  =  Mg(o)  where 


e  6  =  (6/2)Me> 


M  =  £n( 0/2)  +  £n£n(e/2)  +  o(l). 
e 


[Note  that  this  result  only  makes  sense  if  e  -*■  +«>,  i.e.  for  T  <  T  ,  as 

00  a 


supposed. ]  Near  the  turning  point  we  write 


to  obtain 


M  =  M  +  m  with  m  =  0(l) 
e 


D  =  +  2M  (a  +  e"3a)  +  0(l). 


The  function  m  +  e  m  is  graphed  in  Fig.  3.  The  results  (6l)  and  (63),  with 
m  =  0,  are  explicit  asymptotic  formulas  for  the  location  of  the  extinction 
point . 


6.  Other  Aspects  of  Eesnonses. 


We  shall  now  briefly  describe  the  middle  branch  of  the  S-recponse,  the 
two  types  of  monotonic  response  (depending  on  the  nature  of  the  equilibrium 
limit)  and  the  C-response.  A  detailed  treatment  of  all  these  has  been  given 
by  Ludford,  Yannitell  &  Bucknaster  (l976a,b),  in  particular  for  the  C-response 
where  there  is  no  analogous  discussion  of  spherical  diffusion  flames  (Sec.  VI. 6). 


The  flame  temperature  on  the  middle  branch,  which  ranges  from  T  to  T  , 

co  £l 

is  always  the  maximum  temperature  in  the  combustion  field,  so  that  the  reaction 
is  frozen  on  either  side  of  the  flame  sheet.  While  there  could  still  be 
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equilibrium  to  leading  order  beyond  the  flame,  we  shall  first  suppose  that  is 
not  the  case.  As  in  Sec.  VI. 5  we  shall  consider  M  to  be  specified  and 
anticipate  the  equal  conduction  of  heat  to  the  two  sides  of  the  flame.  The 
temperature  profiles  ( VI . 89 )  still  hold  and  lead  to  the  results  ( VI. 90,91) 
with  the  2  deleted.  The  structure  equation  (VI.92)  still  governs,  provided  the 
new  definition 

~e/T# 

(64)  D  =  DLe 

mo* 

is  used.  As  for  the  diffusion  flame,  details  of  the  determination  of  D  are 

m 

not  needed  to  draw  the  conclusion  (VI.95)* 

The  bound  (VI.96)  shows  that  the  solution  holds  for 

(65)  T  <  T  -  1/2 

to  a 

now;  there  are  clearly  no  paramter  values  for  which  it  covers  the  whole  range 
of  M,  since  that  is  infinite  here.  Over  the  rest  6f  the  middle  branch,  i.e. 
the  whole  of  it  when  the  inequality  (65)  is  not  satisfied,  T#  is  zero  and 
the  complete-burning  structure  of  the  monotonic  response  is  needed.  Examples 
of  computed  S-responses  are  shown  in  Fig.  4. 

We  turn  now  to  the  monotonic  responses,  for  which  the  inequality  (VI.99)  holds 
and  the  temperature  increases  beyond  the  flame  sheet.  The  formulas  ( VI. 100, 

101,  102)  are  still  valid,  provided  the  2  is  deleted  from  the  last  two,  so  that 
dTjf/dM  is  again  negative  and  D  is  again  an  increasing  function  of  M  to  leading 
order.  Once  more  the  monotonicity  of  the  response  is  established  without 
discussing  the  flame-sheet  structure.  The  latter  is  of  course  needed  to 
construct  the  curves  themselves  (Fig.  5). 


" - '  — - - 
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As  M  increases  from  its  frozen  value,  r*  decreases  from  00  to  1  and 
there  are  two  possibilities.  If  L  is  less  than  1,  so  that  the  limit  (39) 
does  not  exist,  then  M  increases  without  bound  with  T*  approaching  T&  as 
the  flame  sheet  approaches  the  surface.  The  frozen  limit  is  joined  to*  the 
Burke-Schumann  limit  of  infinitely  rapid  burning.  If  L  is  greater  than  1,  so 
that  the  limit  (39)  does  exist,  then  T#  approaches  Tg  as  the  flame  sheet 
settles  down  on  the  supply  sphere.  The  frozen  limit  is  joined  to  the  Buckmaster 
limit,  the  approach  to  which  is  similar  to  that  for  diffusion  flames  (Sec.  VI. 6). 

Finally  we  come  to  supplies  that  are  hotter  than  the  ambient  atmosphere, 
i  .e.  inequality  (VI.  106).  The  frozen  limit  does  not  exist  but  both  equilibrium 
limits  do.  The  response  takes  the  shape  of  a  C  whose  upper  part  always  corresponds 
to  complete  burning  but  whose  lower  part  nay  correspond  to  incomplete  burning. 

The  division  (if  it  occurs)  lies  below  the  leftmost  point  of  the  C,  which 
corresponds  to  both  ignition  and  extinction  conditions,  so  that  these  conditions 
are  determined  by  the  formulas  (61,63).  An  example  of  such  a  response  is  shown 
in  Fig.  6. 

The  application  of  the  results  in  this  chapter  to  the  burning  of  mono- 
propellant  drops  follows  the  same  lines  as  Sec.  VI. 7.  While  we  have  restricted 
our  attention  to  inert  atmospheres,  the  problem  of  the  ambient  being  an  oxidizing 
atmosphere  for  the  product  of  the  decomposition  should  also  be  investigated.  There 
is  now  a  diffusion  flame  in  addition  to  the  preraixed  one,  so  that  features  from 
Ch.  VI  will  also  arise.  The  Damkohler-number  asymptotics  have  already  been 
considered  by  Fendoll  (19^9)  and,  more  completely,  by  Buckmaster,  Kapila  & 

Ludford  (1978),  but  much  remains  to  be  done. 
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